Symmetry breaking between statistically equivalent, independent channels in a 

few-channel chaotic scattering 
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We study the distribution function P(cj) of the random variable cj = ri/(ri + . . . + tjv), where r k 's 
are the partial Wigner delay times for chaotic scattering in a disordered system with N independent, 
statistically equivalent channels. In this case, r/t's are i.i.d. random variables with a distribution 
>P(t) characterized by a "fat" power-law intermediate tail ~ l/r 1+M , truncated by an exponential 
(or a log-normal) function of r. For N = 2 and N — 3, we observe a surprisingly rich behavior of 
P(u:) revealing a breakdown of the symmetry between identical independent channels. For N — 2, 
numerical simulations of the quasi one-dimensional Anderson model confirm our findings. 
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Scattering by chaotic or disordered systems is encoun- 
tered in various situations ranging from nuclear, atomic 
or molecular physics to mesoscopic devices The key 
property characterizing the scattering process is the uni- 
tary S matrix relating the amplitudes of waves incoming 
to the system and those of scattered waves. Since the 
underlying dynamics is chaotic, the properties of the S 
matrix behave in an irregular way when the parameters 
of the incoming waves or of the medium are modified. 
Hence, an adequate description of the scattering process 
requires the knowledge of the S matrix distribution. 

Time-dependent aspects of the scattering process are 
well captured by the Wigner delay time (WDT) r Q , de- 
fined through the derivative of the S matrix with respect 
to energy E. Physically, r is an excess time spent in the 
interaction region by a wave packet with energy peaked 
at E, as compared to a free wave packet propagation. 

For systems coupled to the outside world via N open 
channels, r = J^feLi r fc; where = dOk/dE are the par- 
tial delay times and 9k are the phase shifts of the S ma- 
trix. One shows that Tfe's are the diagonal elements of 
the Wigner-Smith time delay matrix (WSM), taken in 
the eigenbasis of the scattering matrix 0] . Note that it 
is also customary to define the eigenvalues of the WSM 
as the proper delay times tL (see, e.g., Ref. [HQ). Be- 
yond the 1-channel case, and r' k differ, although their 
sums over all scattering channels are equal to each other. 

Likewise the S'-matrix, the WDT is a random variable, 
whose distribution ^(r) has a generic form pj. IB- 14 1: 
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where a is a characteristic parameter, T(/x) is the gamma 
function and /i is a model-dependent exponent: one en- 
counters situations with < fi < 1, /i = 1 and fx > 1. 

For ID single-channel systems with weak disorder /i = 



1 043, which holds also for quasi- ID disordered sys- 
tems of length L 3> A, where A is the localization length 
Q . One can demonstrate the validity of this result for a 
single-channel scattering in any dimension in the regime 
of strong localization |9|. In ID quasi-periodic systems 
with a single open channel and fractal dimension 
(< 0.5) of the spectrum one has fi = 1 — Dq < 1 [10(, 
and n = 1/2 holds for the 2D generalization of a kicked 



rotor model [11|, [12j , as well as for generic weakly open 
chaotic s yst ems in a parametrically large range of delay 
times [1 IH. Lastly, pL = 1 + Nfi/2 > 1, where f3 is the 
Dyson symmetry index, was obtained for ballistic scat- 
tering from a cavity 0, [TU [lj| ■ 

It is however clear that Eq. (JTJ) defines a limiting form, 
valid either for L — > oo or for weakly open systems. In 
reality, the power-law tail is truncated, such that all mo- 
ments of ^(t) exist. Two mo del- dependent cut-offs seem 
to be physically plausible (although not exact) [H| : 
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K/j,(x) being the modified Bessel function, and a log- 
normally truncated (LNT) form with exp (— In 2 (t)/c) in 
place of exp(— r/b), where b and c are either ~ ^Jl- 
or to the opening degree for weakly open systems |l|, 

In this paper we are concerned with a somewhat un- 
usual statistics of partial delay times for scattering in 
systems with N equivalent channels. We focus here on 
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a random variable which probes the contribution of one 
of the channels to the WDT and hence, the symmetry 
between different channels. To highlight the effect of the 
intermediate power- law tail of \&(r), we suppose that the 
channels are independent of each other such that the par- 
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tial delay times Tfe's are i.i.d. random variables with a 
common distribution in Eqs. ([T]) or © (or a LNT form). 
This situation can be realized experimentally, e.g., for 
scattering in a bunch of disordered fibers. Such a simpli- 
fied model with fx = 1 is also appropriate for a multichan- 
nel scattering from a piece of strongly disordered media 
when the distance between the scattering channels loca- 
tions exceeds A. The role of correlations will be briefly 
discussed at the end of this paper. 

We show here, on example of 2- and 3-channcl sys- 
tems, that intermediate power-law tails entail a surpris- 
ingly rich behavior of the distribution 



P(u) = {S(uj- ri/(n + r 2 + ... + t n ))) 



(4) 



where (...) denotes an average over the distributions of 
Tfe's. We realize that uj exhibits significant sample-to- 
sample fluctuations and, in general, the symmetry be- 
tween identical independent channels is broken, despite 
the fact that all the moments of \&(r) are well defined. A 
similar result was found for related mathematical objects 
in Refs. 15, lal. We address the reader to Ref. Ildl for 



the details on the derivation of P(uj). 

For N — 2 and *( T ) as in Eq. (UJ), we get: 

P{u) = Blu^- 1 (l-ujf 1 . 



(5) 



with B = r(2^)/r 2 (/x). A striking feature of the beta- 
distribution in Eq. ((5]) is that its very shape depends on 
whether < fx < 1, fx = 1 or fx > 1 (see Fig. []}. For 
< (J, < 1, P{ui) is bimodal with a (7-like shape, and 
most probable values being and 1. In this case, the 
symmetry between two identical independent channels is 
broken and either of the two channels provides a dom- 
inant contribution to the WDT. Strikingly, (uj) = 1/2 
corresponds here to the least probable value of uj. For 
fx = 1, P(ui) = 1, and either of the channels may pro- 
vide any contribution to the overall delay time with equal 




FIG. 1. (Color online) N = 2. P(u>) in Eq. © for different 
values of fx. The dashed line depicts P(u>), Eq. ([9]), with /3 = 2 
(,x = 4). 



probability. Finally, for fx > 1, P(uj) is unimodal, which 
signifies that both channels contribute proportionally. 
For N = 2 and a truncated M/(t) as in Eq. @ we get 

P(w) = B[w(l - uj)}- 1 ^ (2y/a/bu>(l-u}j) . (6) 
[2K*J2, 



where B = [2A"^(2 x /a/6)]~ 1 . Note that P{ui) vanishes 
at the edges and is symmetric around uj = 1/2. The 
behavior of P(uj) can be analyzed by expanding the ex- 
pression in Eq. ([6]) in Taylor series at cj = 1/2 [lfj. 
For fx > 1, P(to) is a bell-shaped function with a max- 
imum at uj — 1/2. For fx = 1, for which we previously 
found a uniform distribution, the latter (apart from an 
exponential cut-off at the edges) is approached in the 
limit b/a > 1 [sec Fig[2]a)]. For < fx < 1 the situa- 
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FIG. 2. (Color online) P(u) for N = 2 and *(r) in Eq. © 
with o = 1. a): for jx = 1 and different 6. b): for /x — 1/2 
and different 6. 

tion is more complicated: there exists a critical value b c 
separating two different regimes. For b < b c , P(ui) is a 
bell-shaped function with a maximum at ui = 1/2. For 
b = b c , P(u>) ~ 1 except for narrow regions at the edges, 
where it vanishes exponentially. Finally, for b > b c P(ui) 
has an M-like shape, with maxima close to u> = and 
w = 1, = 1/2 being the least probable value. Hence, 
an exponential truncation of \&(t) does not restore the 
symmetry between different channels, which holds only 
for systems whose size is less than some critical length 
set by b c . Note that for a LNT form the overall behavior 
of P(uj) is very similar and also exhibits a transition at 
some value of the parameter c. 

Further on, for N = 3 and 4 r (r) as in (fTJ), 
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where C 



r(2 M )r(3 M ) 



r »( M )r(2 > i+i/2) and 2*1 is a hyperge- 
ometric series. One finds from Eq. ([7]) that P(uj) ~ 
(1 — cj)' 1_1 when w — > 1 and P(w) ~ w^ -1 when ui — > 0, 
which agrees with the result in Eq. ([5]). On the other 
hand, the amplitudes in these asymptotic forms appear 
to bc very different such that P(uS) is skewed to the left 
[see Fig. [3] a)]. Therefore, for N = 3, P{uj) diverges at 
both edges and has a [/-shaped form for < fx < 1 (sec 
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FIG. 3. (Color online) N = 3. a): P(w) in Eq. © for 
different values of /j,. b): P(uj) in Eq. ((SJ for fi = 1/2 and for 
different values of & (a is set equal to 1). 

Fig. [3] a) , which signifies that the symmetry between the 
channels is broken. For /i > 1 the distribution is uni- 
modal. Remarkably, for /i > 1 the maximum u> m of P(lj) 
is not at oj = 1/3: for [i = 3, one has w m = 0.2719, for 
/i = 10 the maximum is at 0J m = 0.3102 and etc; actually, 
w m — > 1/3 only when fi — > oo. This means that even for 
// > 1, u; exhibits sample to sample fluctuations and the 
average value (oj) = 1/3 docs not have any significance. 
For N = 3 and a truncated ^(t) as in Eq. @ we get 




where J^x) is the Bessel function. We observe that P(oj) 
in Eq. ([5]) is always a bell-shaped function for fi > 1. The 
most probable value oj m is, however, always substantially 
less than 1/3, approaching this value only when fi — > oo 
or b — > 0. The case < /i < 1 is different: For b/a 1, 
P(oj) is peaked at ui m « 1/3. For larger b/a, oj m moves 
towards the origin and P(oj m ) decreases. For yet larger 
b/a, uj m keeps moving towards the origin but now P(oj m ) 
passes through a minimum and then starts to grow. At 
some special value of b/a (b/a « 140 for /i = 1/2) a sec- 
ond extremum emerges (at oj k 0.84 for fi = 1/2) which 
then splits into a minimum and a maximum and P(oj) be- 
comes bimodal. For still larger b/a, the minimum moves 
towards oj = 1/2, while the second maximum moves to 
u = 1. For N = 3 and a LNT form of ^(r), we observe 
essentially the same behavior. 

To substantiate our theoretical predictions, we per- 
formed a numerical analysis of ^(t) and P(u>) for a quasi 
ID disordered Anderson model defined on a rectangular 
lattice of size LxW (with L = 100 and W = 3) with two 
single-channel leads connected to sites (1,2) and [L,2). 
Our (isolated) system is described by the Hamiltonian 
H = X« e »l*)(*l + I*) 01) where Uj are the hop- 

ping rates between the neighboring sites i and j , and ej is 
the energy at the site i, which is a centered, J-correlated 



Gaussian random variable. 

Our numerical results are summarized in Fig. 01 In 
the left panel we depict \1>(t) for different values of the 
localization length A <~ l/(e|). For A/L <C 1, one ob- 
serves that ^(t) decays asymptotically as 1/r 2 , which 
corresponds to fi = 1. On the other hand, ^(r) clearly 
exhibits an intermediate regime with a slower than 1/r 2 
decay (fi < 1). When A increases, this intermediate 
regime shrinks and also the asymptotic decay becomes 
faster (possibly, a log- normal). The right panel shows 
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FIG. 4. (Color online) Distributions vE'(t) and P{w) for a 
quasi-lD Anderson wire for different values of A: from top 
to bottom A = 3, 100 and 200. In the left panel we plot 
r~ 2 (dashed green curve) and r _1 (solid blue curve). In the 
right panel the dashed red curve is the corresponding P(oj) 
calculated under the assumption of statistical independence 
of the TfcS. 

the corresponding distributions P(oj), evidencing a tran- 
sition from J7-like to bell-shaped curves upon an increase 
of disorder. The [/-like shape (A/L <C 1, top right) stems 
out of the intermediate regime with /i < 1 . Interestingly, 
the critical distribution P(oj) ~ 1 is observed for X/L ~ 1 
(middle right), i.e., when A is equal to the length of the 
system. For X/L > 1, a faster than 1/r 2 decay of ^(r) 
leads to a bell-shaped P(w). 

As a test of statistical independence of the actual Tfe's, 
we have computed the distribution P U ncor{^>) (dashed red 
curves in Fig. [J) of the random variable T\/(t\+T2) where 
Ti and T2 are i.i.d. random variables drawn from the nu- 
merically observed ^(t). One notices a good agreement 
between P(lj) (black histogram) and P U ncor{u) , which is 
a clear indication of the lack of correlations between the 
different channels for X/L -c 1. Correlations between 
channels induce some discrepancies between P(to) and 
Puncor(&) only for X/L > 1, when the extension of the 
typical eigenfunction becomes of the order of the system 
size. Consequently, the scattering exhibits a transition 
as the strength of the disorder is varied: T\ and Ti are 
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most likely very different for \/L < 1 and most likely the 
same for X/L > 1. We conjecture that our findings can 
be extrapolated to thin 3D disordered wires, leading to 
a disproportionate contribution of the open channels to 
the total scattering in the diffusive regime, and a propor- 
tionate contribution in the metallic regime. 
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FIG. 5. (Color online) P(u) for the proper delay times for 
different values of /?. a): N = 2, Eq. (JTUJ) . b): N = 3. 

To summarize, we have studied the distribution P(ui) 
of the random variable w, Eq. which defines the con- 
tribution of a given channel to the WDT in a system with 
a few open, independent, statistically equivalent chan- 
nels. We have shown that for 2-channcl systems inter- 
mediate power-law tails with fj, < 1 in the distribution of 
the partial delay times entail breaking of the symmetry 
between the channels; P{uS) has a characteristic £/-shape 
form and the average (ui) = 1/2 corresponds to the least 
probable value. For fi > I the symmetry is statistically 
preserved and (uj) = 1/2 is also the most probable value. 
For = 3 the symmetry between the channels is always 
broken which results in unusual bimodal forms of P(uj). 

Finally, we briefly comment on the effect of correla- 
tions on P{co). We mention two known results on the 
joint distributions of the partial and of the proper delay 
times for which we can evaluate P(w) exactly. The joint 
distribution function of any two partial delay times in a 
system with N channels and arbitrary f3 has been calcu- 
lated in Ref. From this result, we compute exactly 
the distribution P{ui) for two statistically equivalent (but 
not independent) channels: 



P(u) =Duj^I 2 {1-lo) 



(9) 



with D = y*(i+30/2) ' wmcn i s a l so a beta-distribution, 
but with an exponent (= 3/3/2) larger than the one 
(= (3/2) in Eq. (|S|) corresponding to two independent 
channels. For the same f3, the distribution in Eq. ((9]) 
is narrower than P(w) in Eq. (|5|) with /i = 1 + (3/2 
(see Fig. [TJ). Hence, one may argue that the partial de- 
lay times attract each other, which interaction competes 
with the symmetry breaking produced by the intermedi- 
ate power-law tails. Note, as well, that the larger (3 is, 
the narrower is the distribution P(uj). 

The joint distribution of N proper delay times in 



a system with N open channels is also known ex- 
actly Q. It turns out to be given by the Laguerre 
ensemble of random-matrix theory and is defined as a 
product of Ylk=i *( T fc)> where each vE'(r^) as in ([T]) with 
(i = 1 + N/3/2, times the Dyson's circular ensemble, 

Yli<j \ X l T i - 1 / T ']\ f3 - Due t0 the latter factor, the r' k s 
harshly repel each other. For N = 2, we obtain 



P(u J )=Fuj l3 (l-ujy :! \l-2u J \ f3 , 



(10) 



where F is a computable normalization constant. Re- 
markably, P(uj) in ([TO]) is a product of the beta- 
distribution in Eq. ([5]) and a factor |1 — 2uj\p , which is 
a new feature here and stems from the correlations be- 
tween t[ and t-2- This factor forbids r{ and t' 2 to have 
the same values and enhances the symmetry breaking 
[see Fig. flSJ]. Note, however, that two peaks in P{uj) 
become narrower the larger f3 is. Finally, for N = 3, for 
which we can also compute P(ui) exactly, one shows that 
a combined effect of the repulsion and of the intermedi- 
ate power-law tail results in a very peculiar asymmetric 
structure of the distribution [see Fig. ([5])] , which becomes 
increasingly more complicated when (3 increases. 

We thank Y. V. Fyodorov, J. A. Mendez-Bcrmudez 
and T. Kottos for very helpful discussions and comments. 
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